\(\int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx\) [102]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (warning: unable to verify)
   Maple [F]
   Fricas [F(-1)]
   Sympy [F(-1)]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 21, antiderivative size = 78 \[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=-\frac {2\ 2^{5/6} a \operatorname {AppellF1}\left (\frac {1}{2},1,-\frac {5}{6},\frac {3}{2},1-\sin (c+d x),\frac {1}{2} (1-\sin (c+d x))\right ) \cos (c+d x) \sqrt [3]{a+a \sin (c+d x)}}{d (1+\sin (c+d x))^{5/6}} \]

[Out]

-2*2^(5/6)*a*AppellF1(1/2,1,-5/6,3/2,1-sin(d*x+c),1/2-1/2*sin(d*x+c))*cos(d*x+c)*(a+a*sin(d*x+c))^(1/3)/d/(1+s
in(d*x+c))^(5/6)

Rubi [A] (verified)

Time = 0.09 (sec) , antiderivative size = 78, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.190, Rules used = {2866, 2864, 129, 440} \[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=-\frac {2\ 2^{5/6} a \cos (c+d x) \sqrt [3]{a \sin (c+d x)+a} \operatorname {AppellF1}\left (\frac {1}{2},1,-\frac {5}{6},\frac {3}{2},1-\sin (c+d x),\frac {1}{2} (1-\sin (c+d x))\right )}{d (\sin (c+d x)+1)^{5/6}} \]

[In]

Int[Csc[c + d*x]*(a + a*Sin[c + d*x])^(4/3),x]

[Out]

(-2*2^(5/6)*a*AppellF1[1/2, 1, -5/6, 3/2, 1 - Sin[c + d*x], (1 - Sin[c + d*x])/2]*Cos[c + d*x]*(a + a*Sin[c +
d*x])^(1/3))/(d*(1 + Sin[c + d*x])^(5/6))

Rule 129

Int[((e_.)*(x_))^(p_)*((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> With[{k = Denominator[p]
}, Dist[k/e, Subst[Int[x^(k*(p + 1) - 1)*(a + b*(x^k/e))^m*(c + d*(x^k/e))^n, x], x, (e*x)^(1/k)], x]] /; Free
Q[{a, b, c, d, e, m, n}, x] && NeQ[b*c - a*d, 0] && FractionQ[p] && IntegerQ[m]

Rule 440

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[a^p*c^q*x*AppellF1[1/n, -p,
 -q, 1 + 1/n, (-b)*(x^n/a), (-d)*(x^n/c)], x] /; FreeQ[{a, b, c, d, n, p, q}, x] && NeQ[b*c - a*d, 0] && NeQ[n
, -1] && (IntegerQ[p] || GtQ[a, 0]) && (IntegerQ[q] || GtQ[c, 0])

Rule 2864

Int[((d_.)*sin[(e_.) + (f_.)*(x_)])^(n_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Dist[(-b)*(
d/b)^n*(Cos[e + f*x]/(f*Sqrt[a + b*Sin[e + f*x]]*Sqrt[a - b*Sin[e + f*x]])), Subst[Int[(a - x)^n*((2*a - x)^(m
 - 1/2)/Sqrt[x]), x], x, a - b*Sin[e + f*x]], x] /; FreeQ[{a, b, d, e, f, m, n}, x] && EqQ[a^2 - b^2, 0] &&  !
IntegerQ[m] && GtQ[a, 0] && GtQ[d/b, 0]

Rule 2866

Int[((d_.)*sin[(e_.) + (f_.)*(x_)])^(n_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Dist[a^Int
Part[m]*((a + b*Sin[e + f*x])^FracPart[m]/(1 + (b/a)*Sin[e + f*x])^FracPart[m]), Int[(1 + (b/a)*Sin[e + f*x])^
m*(d*Sin[e + f*x])^n, x], x] /; FreeQ[{a, b, d, e, f, m, n}, x] && EqQ[a^2 - b^2, 0] &&  !IntegerQ[m] &&  !GtQ
[a, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\left (a \sqrt [3]{a+a \sin (c+d x)}\right ) \int \csc (c+d x) (1+\sin (c+d x))^{4/3} \, dx}{\sqrt [3]{1+\sin (c+d x)}} \\ & = -\frac {\left (a \cos (c+d x) \sqrt [3]{a+a \sin (c+d x)}\right ) \text {Subst}\left (\int \frac {(2-x)^{5/6}}{(1-x) \sqrt {x}} \, dx,x,1-\sin (c+d x)\right )}{d \sqrt {1-\sin (c+d x)} (1+\sin (c+d x))^{5/6}} \\ & = -\frac {\left (2 a \cos (c+d x) \sqrt [3]{a+a \sin (c+d x)}\right ) \text {Subst}\left (\int \frac {\left (2-x^2\right )^{5/6}}{1-x^2} \, dx,x,\sqrt {1-\sin (c+d x)}\right )}{d \sqrt {1-\sin (c+d x)} (1+\sin (c+d x))^{5/6}} \\ & = -\frac {2\ 2^{5/6} a \operatorname {AppellF1}\left (\frac {1}{2},1,-\frac {5}{6},\frac {3}{2},1-\sin (c+d x),\frac {1}{2} (1-\sin (c+d x))\right ) \cos (c+d x) \sqrt [3]{a+a \sin (c+d x)}}{d (1+\sin (c+d x))^{5/6}} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains complex when optimal does not.

Time = 25.00 (sec) , antiderivative size = 2791, normalized size of antiderivative = 35.78 \[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=\text {Result too large to show} \]

[In]

Integrate[Csc[c + d*x]*(a + a*Sin[c + d*x])^(4/3),x]

[Out]

(3*(a*(1 + Sin[c + d*x]))^(4/3))/(d*(Cos[(c + d*x)/2] + Sin[(c + d*x)/2])^2) - ((15 + 15*I)*AppellF1[2/3, 1/3,
 1/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*(a*(1 + Sin[c + d*x]))^(4/3
)*(1 + Tan[(c + d*x)/2]))/(d*((5 + 5*I)*AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2
- I/2)*(1 + Cot[(c + d*x)/2])]*Sec[(c + d*x)/2] + AppellF1[5/3, 1/3, 4/3, 8/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/
2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*(Csc[(c + d*x)/2] + Sec[(c + d*x)/2]) + I*AppellF1[5/3, 4/3, 1/3, 8/3
, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*(Csc[(c + d*x)/2] + Sec[(c + d*x)/2]
))*(Cos[(c + d*x)/2] + Sin[(c + d*x)/2])^3) + ((15/2 + (15*I)/2)*AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 +
 Tan[(c + d*x)/2]), (1/2 - I/2)*(1 + Tan[(c + d*x)/2])]*(a*(1 + Sin[c + d*x]))^(4/3))/(d*(Cos[(c + d*x)/2] + S
in[(c + d*x)/2])^2*((5 + 5*I)*AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 + Tan[(c + d*x)/2]), (1/2 - I/2)*(1
+ Tan[(c + d*x)/2])] + (AppellF1[5/3, 1/3, 4/3, 8/3, (1/2 + I/2)*(1 + Tan[(c + d*x)/2]), (1/2 - I/2)*(1 + Tan[
(c + d*x)/2])] + I*AppellF1[5/3, 4/3, 1/3, 8/3, (1/2 + I/2)*(1 + Tan[(c + d*x)/2]), (1/2 - I/2)*(1 + Tan[(c +
d*x)/2])])*(1 + Tan[(c + d*x)/2]))) - (3*Cos[(3*(c + d*x))/2]*Csc[c + d*x]*(a*(1 + Sin[c + d*x]))^(4/3)*((1 +
Tan[(c + d*x)/2])/Sqrt[Sec[(c + d*x)/2]^2])^(2/3)*(8 + (1 + I)*2^(2/3)*(((1 - I)*(I + Cot[(c + d*x)/2]))/(1 +
Cot[(c + d*x)/2]))^(1/3)*Hypergeometric2F1[1/3, 2/3, 5/3, ((1 + I) + (1 - I)*Tan[(c + d*x)/2])/(2 + 2*Tan[(c +
 d*x)/2])]*(I + Tan[(c + d*x)/2]) - AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/
2)*(1 + Cot[(c + d*x)/2])]*((2 + 2*I) - (2 - 2*I)*Cot[(c + d*x)/2])^(1/3)*((-1 - I)*(I + Cot[(c + d*x)/2]))^(1
/3)*(1 + Tan[(c + d*x)/2])))/(4*d*(Cos[(c + d*x)/2] + Sin[(c + d*x)/2])^3*(1 + Tan[(c + d*x)/2])*((-3*Sec[(c +
 d*x)/2]^2*((1 + Tan[(c + d*x)/2])/Sqrt[Sec[(c + d*x)/2]^2])^(2/3)*(8 + (1 + I)*2^(2/3)*(((1 - I)*(I + Cot[(c
+ d*x)/2]))/(1 + Cot[(c + d*x)/2]))^(1/3)*Hypergeometric2F1[1/3, 2/3, 5/3, ((1 + I) + (1 - I)*Tan[(c + d*x)/2]
)/(2 + 2*Tan[(c + d*x)/2])]*(I + Tan[(c + d*x)/2]) - AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*
x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*((2 + 2*I) - (2 - 2*I)*Cot[(c + d*x)/2])^(1/3)*((-1 - I)*(I + Cot[
(c + d*x)/2]))^(1/3)*(1 + Tan[(c + d*x)/2])))/(8*(1 + Tan[(c + d*x)/2])^2) + ((8 + (1 + I)*2^(2/3)*(((1 - I)*(
I + Cot[(c + d*x)/2]))/(1 + Cot[(c + d*x)/2]))^(1/3)*Hypergeometric2F1[1/3, 2/3, 5/3, ((1 + I) + (1 - I)*Tan[(
c + d*x)/2])/(2 + 2*Tan[(c + d*x)/2])]*(I + Tan[(c + d*x)/2]) - AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 +
Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*((2 + 2*I) - (2 - 2*I)*Cot[(c + d*x)/2])^(1/3)*((-1 - I
)*(I + Cot[(c + d*x)/2]))^(1/3)*(1 + Tan[(c + d*x)/2]))*(Sqrt[Sec[(c + d*x)/2]^2]/2 - (Tan[(c + d*x)/2]*(1 + T
an[(c + d*x)/2]))/(2*Sqrt[Sec[(c + d*x)/2]^2])))/(2*(1 + Tan[(c + d*x)/2])*((1 + Tan[(c + d*x)/2])/Sqrt[Sec[(c
 + d*x)/2]^2])^(1/3)) + (3*((1 + Tan[(c + d*x)/2])/Sqrt[Sec[(c + d*x)/2]^2])^(2/3)*(-1/2*(AppellF1[2/3, 1/3, 1
/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*((2 + 2*I) - (2 - 2*I)*Cot[(c
 + d*x)/2])^(1/3)*((-1 - I)*(I + Cot[(c + d*x)/2]))^(1/3)*Sec[(c + d*x)/2]^2) + ((1 + I)*(((1 - I)*(I + Cot[(c
 + d*x)/2]))/(1 + Cot[(c + d*x)/2]))^(1/3)*Hypergeometric2F1[1/3, 2/3, 5/3, ((1 + I) + (1 - I)*Tan[(c + d*x)/2
])/(2 + 2*Tan[(c + d*x)/2])]*Sec[(c + d*x)/2]^2)/2^(1/3) + ((1/3 + I/3)*2^(2/3)*(((1/2 - I/2)*(I + Cot[(c + d*
x)/2])*Csc[(c + d*x)/2]^2)/(1 + Cot[(c + d*x)/2])^2 - ((1/2 - I/2)*Csc[(c + d*x)/2]^2)/(1 + Cot[(c + d*x)/2]))
*Hypergeometric2F1[1/3, 2/3, 5/3, ((1 + I) + (1 - I)*Tan[(c + d*x)/2])/(2 + 2*Tan[(c + d*x)/2])]*(I + Tan[(c +
 d*x)/2]))/(((1 - I)*(I + Cot[(c + d*x)/2]))/(1 + Cot[(c + d*x)/2]))^(2/3) - ((1/6 + I/6)*AppellF1[2/3, 1/3, 1
/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*((2 + 2*I) - (2 - 2*I)*Cot[(c
 + d*x)/2])^(1/3)*Csc[(c + d*x)/2]^2*(1 + Tan[(c + d*x)/2]))/((-1 - I)*(I + Cot[(c + d*x)/2]))^(2/3) - ((1/3 -
 I/3)*AppellF1[2/3, 1/3, 1/3, 5/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2])]*((-
1 - I)*(I + Cot[(c + d*x)/2]))^(1/3)*Csc[(c + d*x)/2]^2*(1 + Tan[(c + d*x)/2]))/((2 + 2*I) - (2 - 2*I)*Cot[(c
+ d*x)/2])^(2/3) - ((2 + 2*I) - (2 - 2*I)*Cot[(c + d*x)/2])^(1/3)*((-1 - I)*(I + Cot[(c + d*x)/2]))^(1/3)*((-1
/30 + I/30)*AppellF1[5/3, 1/3, 4/3, 8/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 - I/2)*(1 + Cot[(c + d*x)/2]
)]*Csc[(c + d*x)/2]^2 - (1/30 + I/30)*AppellF1[5/3, 4/3, 1/3, 8/3, (1/2 + I/2)*(1 + Cot[(c + d*x)/2]), (1/2 -
I/2)*(1 + Cot[(c + d*x)/2])]*Csc[(c + d*x)/2]^2)*(1 + Tan[(c + d*x)/2]) + ((2/3 + (2*I)/3)*2^(2/3)*(((1 - I)*(
I + Cot[(c + d*x)/2]))/(1 + Cot[(c + d*x)/2]))^(1/3)*(I + Tan[(c + d*x)/2])*(2 + 2*Tan[(c + d*x)/2])*(-((Sec[(
c + d*x)/2]^2*((1 + I) + (1 - I)*Tan[(c + d*x)/2]))/(2 + 2*Tan[(c + d*x)/2])^2) + ((1/2 - I/2)*Sec[(c + d*x)/2
]^2)/(2 + 2*Tan[(c + d*x)/2]))*(-Hypergeometric2F1[1/3, 2/3, 5/3, ((1 + I) + (1 - I)*Tan[(c + d*x)/2])/(2 + 2*
Tan[(c + d*x)/2])] + (1 - ((1 + I) + (1 - I)*Tan[(c + d*x)/2])/(2 + 2*Tan[(c + d*x)/2]))^(-1/3)))/((1 + I) + (
1 - I)*Tan[(c + d*x)/2])))/(4*(1 + Tan[(c + d*x)/2]))))

Maple [F]

\[\int \csc \left (d x +c \right ) \left (a +a \sin \left (d x +c \right )\right )^{\frac {4}{3}}d x\]

[In]

int(csc(d*x+c)*(a+a*sin(d*x+c))^(4/3),x)

[Out]

int(csc(d*x+c)*(a+a*sin(d*x+c))^(4/3),x)

Fricas [F(-1)]

Timed out. \[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=\text {Timed out} \]

[In]

integrate(csc(d*x+c)*(a+a*sin(d*x+c))^(4/3),x, algorithm="fricas")

[Out]

Timed out

Sympy [F(-1)]

Timed out. \[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=\text {Timed out} \]

[In]

integrate(csc(d*x+c)*(a+a*sin(d*x+c))**(4/3),x)

[Out]

Timed out

Maxima [F]

\[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=\int { {\left (a \sin \left (d x + c\right ) + a\right )}^{\frac {4}{3}} \csc \left (d x + c\right ) \,d x } \]

[In]

integrate(csc(d*x+c)*(a+a*sin(d*x+c))^(4/3),x, algorithm="maxima")

[Out]

integrate((a*sin(d*x + c) + a)^(4/3)*csc(d*x + c), x)

Giac [F]

\[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=\int { {\left (a \sin \left (d x + c\right ) + a\right )}^{\frac {4}{3}} \csc \left (d x + c\right ) \,d x } \]

[In]

integrate(csc(d*x+c)*(a+a*sin(d*x+c))^(4/3),x, algorithm="giac")

[Out]

integrate((a*sin(d*x + c) + a)^(4/3)*csc(d*x + c), x)

Mupad [F(-1)]

Timed out. \[ \int \csc (c+d x) (a+a \sin (c+d x))^{4/3} \, dx=\int \frac {{\left (a+a\,\sin \left (c+d\,x\right )\right )}^{4/3}}{\sin \left (c+d\,x\right )} \,d x \]

[In]

int((a + a*sin(c + d*x))^(4/3)/sin(c + d*x),x)

[Out]

int((a + a*sin(c + d*x))^(4/3)/sin(c + d*x), x)